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Abstract 

We have studied a DGP-inspired braneworld scenario where the idea of Lorentz 
invariance violation has been combined into a specifying preferred frame that em- 
bed a dynamical normal vector field to brane. We propose the Lorentz violating 
DGP brane models with enough parameters can explain crossing of phantom divide 
line. Also we have considered the model for proper cosmological evolution that is 
according to the observed behavior of the equation of state. In other view point, 
we have described a Rip singularity solution of model that occur in this model. 
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1 Introduction 



According observation data from WMAP, universe accelerated expanding at the current 
epoch [1,2,3,4,5]. These astrophysical observations also show that the universe is spatially, 
and is composed of about 27% of dark matter, and 73% of homogeneously distributed 
new type of negative pressure matter, known as dark energy, that lead to the current 
accelerated expansion in universe. From another view point, theories of extra dimensions, 
in which the observed universe is realized as a brane embedded in a higher dimensional 
spacetime, have attracted some attentions. Particularly, the model proposed by Dvali, 
Gabadadze and Porrati (DGP) [6,7,8] is different from others. In this regard, DGP model 
also predicts deviations from the standard 4-dimensional gravity in over large distances. 
However, impact of Lorentz invariance violating models (LIV) on cosmology has been 
studied in [9,10,11]. This models has been used in the context of scalar- vector-tensor 
theories [9,11]. One important observation has been made in reference [11,12] which 
considered accelerated expansion and crossing of phantom divide line with one minimally 
coupled scalar field in the attendance of a Lorentz invariance violating vector field. As 
one has shown in [13], quintessence model with a single minimally coupled scalar field 
has not the capability to describe crossing of phantom divide line, uj — — 1. But, a single 
non- minimally coupled scalar field is suitable to cross the phantom divide line [14]. 

In following, we study a new method of performing local Lorentz violation in a grav- 
itational setup to make the existence of a tensor field with a non-vanishing expectation 
value, and then add this tensor to gravity or matter and scalar fields. The example of 
this approach is to used a single time-like vector field with unit norm. This vector field 
choose a preferred frame at each point in space-time and a matter field coupled to it will 
endure a violation of local Lorentz invariance. A method of this model was introduced 
Kostelecky and Samuel in [15]. Here we choose the vector to have dynamics, and fixed 
norm in the action. 

Our procedure in this paper is deriving the basic equations of motion for the most 
general theory of a time-like vector field w M with an action S = Ssuik + Ssrane where 
S Brane = Seh + S<j) + S m + S u [11]. Then with numerical calculations, we study some 
cosmological aspects in this model in present the vector field in our formalism. In this 
regard, we study other solutions admitted as Rip singularity , that occur in the condition 
u < — 1 increases rapidly. However, it possible different types of singularity, depending 
energy density and scale factor how increases with time[16-19, 31-36]. In following we just 
focus on Big or Little singularity and other cases of singularity maybe appear in this 
model with tuning parameters space. 
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2 A Lorentz Violating DGP Brane Model 



As we shown in[ll], for study Lorentz violation, we assume a vector in the extra 
dimension. So, this local frame at space-time is unavoidably selected as the preferred 
frame. However the action of the Lorentz violating DGP scenario in the presence of a 
minimally coupled scaler field and vector field on the Brane can be written as the sum of 
two separate parts 

S = Ssulk + SBrane, (1) 

where S Bu ik and S Bra ne = S EH + + S m + S u are 

/ jji ^ 
d 5 x^V=9K, (2) 

J d 4 x^(^R[q} - ^<T V^V„0 - V(</>) + mlK + £ m + 



''Brane 



V WVnUy-fo VVV„v^3 (V^^-A/u^V./V.Ua+A + 1)] J . (3) 

7 -I y=o 

y is coordinate of the fifth dimension, m\ point to the constant in the Bulk and m§ used 
for Brane, for more details see [11]. We consider brane is located at y — O.qab is five 
dimensional bulk metric with Ricci scalar 1Z, also q^ v is induced metric on the brane with 
induced Ricci scalar R. g^ B and q^ v are connected via = 5 U & V B Qab- This action is 
suitable for add any non-gravitational degrees of freedom in the model of Lorentz violating 
scalar-tensor- vector theory. In other hand, we take u^u u = — 1 and the expectation value 
of vector field -u M is < 0|-u M -u^|0 >= — 1 [20]. (i = 1, 2, 3, 4) are temporary parameters 

with dimension of mass squared, A is also a Lagrange multiplier field. Note that are 
mass scale of violating Lorentz invariance [9,20,21]. 

A motivation for choosing the vector action in term (3) discussed in Ref. [20]. Also 
cosmological consequences of this action are studied in Ref. [9,11,21]. 
In the action (3) ordinary matter part is shown by Lagrangian C m = Cm^q^^ip) where ip 
is matter field and its energy- momentum tensor is [22] 

Tfj, v = — + q^C-m. (4) 

The scalar field Lagrangian, = — a 4>V u4> ~ ^(0)> gives the following energy- 
momentum tensor 

V = V M 0V,0 - -^„(V0) 2 - quuV{4>). (5) 
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The energy-momentum tensor of vector field also is by usual formulate 

5C^ 



rp(u) _ _9. 

8q" v 



(6) 



The Bulk-brane Einstein's equations obtained from action (1) are given by 

ml (ll AB - l -g AB TZ) + m\5 A »5 B v (r^ - ^i?) 5{y) = 5 A W%) 



(7) 



where is 4-dimensional(brane) d'Alembertian and T^ u — , , ^ 
From equation (7) we find 



T _|_ r _|_ t(u) 



G AB = 1Z AB — -g AB lZ = 



and 



— I R 



T 



(8) 



(9) 
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for bulk and brane respectively. We use the following line element to derive cosmological 
equations of our model 

ds 2 = q^dx^dx" + b 2 (y, t)dy 2 = -n 2 (y, t)dt 2 + a 2 (y, t)% j dx i dx j + b 2 (y, t)dy 2 (10) 

we use n(0,t) = 1 and Af = 1 after the variation [11]. 

Now we can obtain components of the total energy-momentum tensor as follow 



and 



P = Pm + Pu + P<t> 



P=Pm+Pu+P<t> 



(11) 



;i2) 



Here we take the energy-momentum tensor for the matter as a perfect fluid with energy 
density p m and pressure p m is T^ u = (p m + p m )N^N u + p m q^ u , where is a unit time-like 
vector field illustrate the fluid four-velocity. We also assume a usual equation of state for 
the fluid as p m = (7 — l)p m that 1 < 7 < 2. Energy density and pressure of minimally 
coupled scalar field are taken as 



1 • 



(j) 2 + n 2 V((f)) 



and 



P<t> 



2n 2 



2 - Vtf) 



y=0 



(13) 



(14) 
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where a dot point to the derivative respect to cosmic time t. The stress-energy for the 
vector field also assume the form of a perfect fluid, with an energy density determined by 
[10] 

Pu = -3(3H 2 (15) 

and a pressure 



p u = (5H 2 \3 + 2^- + 2 



H 2 H/3 



(16) 



where we have assume the parameter /3 = (ft + 3ft + ft) and H = £2j*g is Hubble 
parameter. In the absence of vector field, all ft = 0, In following, we stress for obtain 
equation of state in general case one should varying the action (1), that It has been calcu- 
lated in Ref[9,10 11,12,15,20, 21 and 22]. Here we just used their results and summarized 
calculations. Hence, we write the general equation of state which takes the following form 



3 2{fr ■ 2- 



+ - n<t>) 



Pm -3fttf 2 + k> 2 + n 2 \/(0) 



2/=0 



(17) 



and for obtain generalized Friedmann equations, we use the effective Einstein equation 
(for more details see[ll,23]), therefore we have 



12m*j 



Pr, 



3(5 H 2 + U 2 + nV(0) - 3m 2 3 (H 2 + 



(18) 



As we calculated in above equations, there are many parameters, so we expect the Lorentz 
violating brane model, suitable for study the current stage acceleration of the universe 
and other cosmological aspects with fine tuning parameters. 



3 Setup for useful equations 



In following, we consider scalar equation of state u^t) which is a candidate for the dark 
energy in future discussions. In this regard, the equation of state parameter on the brane 
given by 



UJ4, = 



P<t> 



^ 2 ~ v(4>) 



U 2 + n 2 V((p) 



(19) 



Now we obtain scalar field of equation of state in two ways. First, we have Friedmann 
equation from equation (18) as 



H 



So 1 

3a 4 36m| 



- 3(3H 2 + \$ + V(<j>) - 3m 2 3 H 2 



(20) 
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where k = and we assume non matter contain on the model that means p m = and 
p m = 0, where £° + AH£ = and £° = ^ with £ as an integration constant [23]. 
Then we obtain dynamic of scalar field from above equation 



2 = 6tf 2 (/3 + ml) - 2V(<f>) + Umle^H* + ^ (21) 

where e = ±1 has main role in next calculations. For determine dynamic of scalar field, 
we use to a well-known potential: models with potential of the form V((f>) = \'<p 2 [24]. 

Second, Such as we see equation (21) depend on potential functions, therefore we will 
obtain the scalar field of equation of state in our model that independent of potentials. 
So the energy equation for vector field taken as 

p u + 3H( Pu + Pu ) = +3H 2 p (22) 

and for the scalar field 

p 4> + 3H{p 4> +p ( , ) ) = -3H 2 $. (23) 
The total energy equation in presence of both the vector and the scalar fields is 

p + 3H(p + p)=0, (p = p u + P4> ) (24) 

Now we obtain dynamic of scalar field by differentiating equation (13) with respect to t 
and using equation (23) then we have 

4> + 3Hj> + 3H 2 (5^ + Vj = (25) 

and by differentiating equation (20) with respect to t and using equation (25) we have 

= T2ml( 2 -§^ + -g^)* - 2Hfij - 2fiHj - (26) 
<pH <p z aa A ti 

where we assume H and j3 depended on . For simplicity in following we take £ = 0. 
In other hand, with substituting Equation (26) into the Friedmann equation (20) obtain 
the potential of the scalar field as 

Tr/J , 1 12H + 6m 3 2 # 2 m 4 3 - 2 m 4 3 + 6H 2 m 4 3 /3 

V{<P) = 2 W (27) 



and 

1 -12H + 6m 3 2 H 2 m 4 3 - 2 m 4 3 + 6^ 2 m 4 3 /3 

2 m 4 



v{4>) = ;"" liTU " t3 - ( 28 ) 
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1 -10 J 



sing= 1, zeta= 1, xi= 1 sing= -1, zeta= 1, xi= 1 

sing= 1, zeta= -1, xi= -1 sing= -1, zeta= -1, xi= - 



Figure 1: Variation of ^ relative to 4> in equation (29) for different values of £ and ( according 
ansatz H = Hq^ and (3(cf)) = m(jfi. 

Also factor from equation (26) give as 

= 1 + 2/9 /9if + .6/3^ + .6/3#x 

(I (lOH 4, + 4 + 40H 4, + 5/3^ + 4/3^/3 + + 4) 

^ : — ^ : (29) 

xC(0) 

where A = 2^100 + 10H^ + 30/3 + 40/3 + 30/3^/3 + 10/3# e' H and e' = ± 1. We em- 
phasize that equation (26) have several solutions, here we choose a real solution of them. 
Besides, equation (29) has not the analytical solution that we exactly use it for cosmo- 
logical consideration. Hence, we just study some phenomenological aspects from Figure 
1 and choose a best fit scalar filed with equation (29) for analytical consideration. 
In follows, we calculate scalar field of equation of state with using equation (20). Therefore 
we have two equations for if 2 as 

2 1 P<p ^3 2 + P P<p + 6 m 4 6 + 2 ^3m 4 % m 3 2 + 3 m 4 6 /3 + 9 m 4 12 
H = 3 2/3m 3 2 + /3 2 + m 3 4 (30) 

and 

2 1 P4> m 3 2 + /3 + 6 m 4 6 - 2 ^3m 4 % m 3 2 + 3 m 4 6 /3 p + 9 m 4 12 
~3 2/3m 3 2 + /3 2 + m 3 4 ' ^ 

For each branch of if 2 means equation (30) and (31) with using equation (23) we obtain 
equation 

El + E2(l +0J4) = -$E3 (32) 



7 



where 



and 



E2 



Ptf, H (3 Hm 4 + (3 H m 4 3 + m 2 H m 4 3 

frj I — — ~|- 

p H 2 e 4 + (3 Hm A 3 + m 3 2 Hm 4 3 



(/? + m 3 2 ) 



(33) 



(34) 



and 



E3 = 



2 e 3 ^m 4 3 (m 4 9 + iJ 2 (/3 + m 3 2 ) 2 m 4 3 + 2 #e 4 (/3 + m 3 2 ))m 4 s 
# (m 4 3 05 + m 3 2 ) + 2e 4 ) (/3 + m 3 2 ) 2 

# 2 (/3 + m 3 2 f m 4 3 + 2 He A (f3 + m 3 2 ) + 2 m 4 9 



(35) 



# (m 4 3 (f3 + m 3 2 ) /J + 2 e 4 ) (/3 + m 3 2 )' 
where e\ = e 2 = e 3 = e 4 = ±1. We can calculate after simplify from equation (32) as 



u<j> = 1/6 



- 6 e 2 (/3 + m 3 2 ) #ei (l/2 m 4 3 (/3 + m 3 2 ) # + e 4 ) x 



1 



m 4 9 + H 2 ((3 + m 3 2 ) 2 m 4 3 + 2 ffe 4 Qg + m 3 2 ) 
m 4 3 



2 e 3 ^/m 4 3 (m 4 9 + H 2 ((3 + m 3 2 ) 2 m 4 3 + 2 #e 4 (/3 + m 3 2 ))m 4 3 - 
3 m 4 3 (/3 + m 3 2 ) 2 (2/3 ^ + e 1 m 4 3 ) # 2 - 



6 ((l/3 ify m 3 4 + 2/3 # /5 m 3 2 + Cl e 4 + 1/3 /3 2 # 0) m 4 3 + 1/3 ^ e 4 ) (/3 + m 3 2 ) 



2 (m 4 9 /3 + #^e 4 (/3 + m 3 2 )' 



m 4 + e 2 



V 



m 4 9 + H 2 (13 + m 3 2 ) 2 m 4 3 + 2 #e 4 (/3 + m 3 2 ) 



m 4 ^ 



x 



(/3 + m 3 2 ) 1 H- X e x - X (l / 2 m 4 3 (f3 + m 3 2 ) H + e 4 ) 



(36) 



In following we study our models by using above equations for obtain some constraint of 
parameters that be according to observation data and also discuss about other cosmolog- 
ical aspects. 
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4 Fine tuning parameters 



In this section we consider our equation that obtained in above section for determine 
some cosmological aspects. Therefore we need to solve equation (36) to study crossing of 
phantom divide barrier in model. In first stage we should obtain dynamics of scalar 
field with equation (21). This will be achieved only if the Hubble parameter H(<f)(t)) 
and the vector field coupling, /3((f>(t)) are known. Hence, our strategy is to choose some 
cases of the Hubble parameter H(<f)(t)) and the vector field coupling /3(<f>(t)) for considering 
possible crossing of phantom divide barrier (PDL) and Rip singularity in this model. Also 
we should obtain suitable domains of parameter space which have the capability to explain 
Rip singularity and crossing of phantom divide line by equation of state parameter. 

In following, we take a general case of the vector field coupling and the Hubble pa- 
rameter that are functions of scalar field 

H = H (f)< , /3(<f>) = mtf (37) 

where Hq and m are positive constant parameters. Here we used anzats (37) because 
considered some cases of the solution and verified the stability in previous works [12]. Also 
other authors [10 and references therein] have been used anzats (37) for some cosmological 
solutions, for example, late time acceleration, deceleration parameters and etc. So let us 
using equation (21) for obtain dynamic of scalar field equation with assuming case (37) 
as 

= 6# O 2 2C (™<^ + m 3 2 ) - X'(f) 2 + 12m 4 3 e# o f (38) 
for potential V(4>) = X'cj) 2 . 

If we use a simple case from a acceptable range of £ and ( for analytical solution of scalar 
field equation (38), we obtain analytical solution for dynamic of scalar filed equation as 
follows 

1 



m 2 



12e ( t+ A )V6tfoW-A' J 6Ho 2 m3 2 _ X 'm 4 3 eH + 6 e »v^oW^ : H 2 m 3 2 



(6/JoW-A') 3/2 ] _1 (39) 



-0+A () )A/6H 2 m 3 2-A' 
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where we take £ = — 1 , £ = 1 and A is an integration constant. 

In the following we consider equation of state (36) with using equation (39) for exact 
dynamic of equation of state, Then study crossing phantom divided line and other options 
in our model. The result of our numerical calculation has been shown in figures 2-5 and 6-9. 
The summarize of results given in table 1. According figures 2-5 and table 1, Our model 
can explain crossing phantom divided, that is relevant to observational data. In other 
hand, one has been discussed in [29-36] about what conditions lead to a Rip singularity 
solutions, now we consider, can this model predict a Rip singularity in context, we note 
that just two cases of singularity means Big and Little singularity study in following. 

Therefore, as shown in figures 6-9, we provide suitable conditions for appear Rip 
singularity (this situation occur in e = +1 branch). As see in a numerical solution, this 
situations appear in the phantom barrier, in this case the equation of state increase rapidly 
in finite or infinite time. Hence we can consider the solution in two major cases as known 
Little and Big singularities. In fact we can take a dynamical explication to a case of 
solutions that leading to the requirement of Big Rip or Little Rip singularity. Instead, 
from this case of solutions we should choose those requirement which the constraint by 
observational data. In other view point, some other types of singularity maybe appears 
in our work, one can easily check other type singularity with suitable conditions in this 
model. 

Now we should describe what case of equation of state is studied for determine the char- 
acter of mixed fluids. Typically in this model, we use three energy-momentum contents, 
1: The ordinary matter 2: Dark energy describe by the scalar field 3: energy-momentum 
content which determine by a Lorentz violation vector field. Here we consider that ordi- 
nary matter has partial part of entire energy- momentum content. But for other contents, 
in our model it is possible to take the "trigger mechanism" to describe dynamical equa- 
tion of state, i.e. we suppose that scalar- vector-tensor theory with LIV acts similar a 
hybrid inflation models. In this regard, scaler and vector field have the roles of inflaton 
and the "waterfall" field. Hence, we can fine-tuning parameters of model for best fit with 
observational data [10,25]. It is acceptable to expect that one of them suddenly dominate 
and we have a cosmological stage, for example, inflation phase or acceleration phase and 
etc. Note that an important result in this context is crossing of phantom divide barrier 
PDL that determine with a single minimal coupling scalar field and by a Lorentz violating 
vector field, if considering suitable fine tuning of model parameters. 
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Figure 2: Crossing of phantom divide line. 
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Figure 3: Crossing the phantom divide line. 
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Figure 4: Crossing the phantom divide line. 
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Figure 5: Crossing the phantom divide line. 
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Figure 7: Rip singularity 
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Figure 8: Rip singularity 



Figure 9: Rip singularity 



12 



Table 1: Considering Crossing phantom divided line and Rip singularity in equation of 
state in equation (36) on all branch for e = — 1 (middle) and e = +1 (right) in equation 
(21) for potential of the form V(4>) = A'0 2 



ei ^2 ^3 ^4 


Crossing PDL 


Rip singularity 


+ + + + 


Yes 


Big 


+ + + - 


Yes 


Big 


+ + - - 


Yes 


Little 


+ - - - 


Yes 


Big 




Yes 


normal 


- + + + 


Yes 


Big 


h + 


Yes 


Big 


1_ 


Yes 


normal 


+ - + + 


Yes 


crossing PDL 


H h 


Yes 


Big 


+ + - + 


Yes 


normal 


+ - + - 


Yes 


Big 


•• + - + 


Yes 


Big 


- + + - 


Yes 


Big 
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5 Summary 



We have consider a possible violation of Lorentz invariance in a DGP Brane cosmology. 
Also we have shown that by a suitable choice of parameter space, it is possible to have 
crossing phantom divided barrier (PDL) and it can predicts Rip singularity in LIV context. 
We used an interactive picture, a minimally coupled scalar field and a Lorentz violating 
vector field can lead to the phantom phase. The comparison of our results with cosmology 
of 4D non-minimal vector theories show same aspects but with non minimal coupling in 
vector term, for example, explain the late time acceleration of the Universe [26,27,28]. In- 
stead, those models have been shown that a non minimally vector field can be a candidate 
for dark energy [28] just like our model. 

In other view point, in this framework, there is the possibility of a Rip singularity by 
suitable tuning in the parameters. As shown in Figs. 6-9, we have studied some type of 
solutions called Little or Big Rip, that maybe occur in structures of the Universe, accord- 
ing to the increasing of the dark energy density. However, In correct choice of parameters 
maybe other types of singularity appear in this model. Finally we emphasize that we 
have used a model of a scalar field with Lorentz violation to describe some cosmological 
aspects. We believe this model can perform some fascinating predictions for observations. 
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